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Techniques in High-Resolution Coincidence Counting 

George H. Minton ' 

Circuitry, systems, and techniques used in radiation coincidence measurements with 
millimicrosecond resolution are discussed. An analysis of the time fluctuations in a scintilla- 
tion counting system is presented and comparisons with experimental results are made. 



1. Introduction 

The applications of coincidence systems to prob- 
lems in nuclear research are of two types. In one, 
the coincidence system is used in order to detect 
events that arc coincident in time. In the other, it 
is applied to the measurement of the distribution in 
time of noncoincident events. The characteristics 
of some coincidence systems have been studied in 
order to determine their suitability for such applica- 
tions. The factors that have been considered are 
(1) the resolution obtainable with available radia- 
tion detectors, (2) the relationship between the reso- 
lution and the efficiency of the system, (3) the fac- 
tors that determine the response of the system to 
coincident events and to events separated in time by 
intervals of the order of magnitude of the resolving 
time or less. 

Although the ideal coincidence system will respond 
only to events that are exactly coincident, such 
response is not obtainable in practice. The events 
involved arc signaled by the production in the de- 
tector of electrical pulses — often widely varying in 
amplitude — and the resolution of the system is de- 
termined by the response of the detector and that of 
the circuit to which the signals are fed. Spurious 
response of such a system corresponds to the re- 
cording as coincidences of single pulses or of events 
that are more closely separated in time than its 
resolution. As the number of the latter events in 
any given time interval depends on the product of 
the counting rates in the individual detectors, the 
resolving time of the system often provides a limita- 
tion on the permissible counting rates. Failure of 
the system to respond to coincident events can be 
measured in terms of an efficiency factor that is less 
than or equal to 1. 

In the measurement of the distribution in time 
of noncoincident events, the signal from one detector 
is delayed before being fed into the coincidence cir- 
cuit, and the number of coincidence's recorded is 
studied as a function of the length of the inserted 
delay time. These data furnish a response curve of 
the system for the events under study. If the rate 
of occurrence of events does not vary appreciably in 
a time of the order of the resolution of the system, 
the shape of the response curve of the system to 
coincident events is not critical. On the other hand, 
for studying more rapidly varying distributions it 
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is necessary to obtain a response curve of the sys- 
tem for coincident events. The resolution of the 
system is often defined as one-half the time as 
measured on the delay axis between points cor- 
responding to one-half the maximum counting rate 
of such a response curve. By proper analysis of 
this curve along with that for noncoincident events, 
it is possible to measure time intervals to an order 
of magnitude smaller than the 4 resolution of the 
system, as defined above, to an accuracy of the 
order of 10 percent. 

2. Detectors 

The only type of radiation det eel or presently 
available that is suitable for very fast coincidence 
work is the scintillation counter. Among the vari- 
ous phosphors available, there is a wide divergence 
among the values of the lime constants for the emis- 
sion of light. In general, the organic phosphors (unit 
their light in one-tenth to one one-hundredth the lime 
of the inorganic materials, so that attention will be 
confined to the scintillation counter with organic 
phosphor. The measurements reported here were 
made, using RCA 5819 and 0199 photomultipliers 
with stilbene crystals or solutions of terphenyl in 
toluene as scintillators. 

3. Coincidence Circuit 

The coincidence circuit receives the pulses from 
the detectors and provides an output pulse whenever 
the input pulses are in coincidence. The ideal co- 
incidence circuit never responds to single pulses, re- 
gardless of their amplitude, and never registers a co- 
incidence when the input pulses are separated by a 
time greater than their own duration. The circuit 
should have as low an input capacitance as possible 
in order to minimize any effect on the rise time of the 
input pulse. The requirement of low capacitance is 
most satisfactorily met by the use of a circuit con- 
structed from crystal diodes. Figure 1 is a schematic 
diagram of such a circuit — a modification of that of 
DeBenedetti and Richings [1] 2 — which has been used 
in this investigation. Study of the response of this 
circuit to pulses from a mercury-relay pulse generator 
and from scintillation counters has revealed that the 
response of the circuit is not a limiting factor in the 
operation of the coincidence systems using nuclear 
radiation detectors, so that an extensive comparison 
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Figure 1. High-resolution coincidence circuit for negative 
pulses from the anodes of two photomultiplier tubes. 




5 6 

INPUT NO. I. 



Figure 2. Response of the fast coincidence circuit of figure 1 and 
associated amplifier to simultaneous pulses from a mercury- 
pulse generator. 



Fast-rising negative pulses having an exponential decay time of 2.5xl0~ 9 sec 
were generated by discharging a capacitor through the relay into a coaxial cable. 
The pulses were further shaped, as in the scintillation counters, by reflection in 
5-in. shorted stubs of RG 62-U coaxial cable. 

of available coincidence circuits has not been made. 
On the basis of published information, however, it is 
believed that this circuit is at least as good as any 
other in use today. Figures 2 and 3 illustrate the 
response of this circuit to pulses from the mercury- 
relay pulse generator. In order to closely approxi- 
mate the shape of the pulse from a scintillation 
counter with liquid scintillator, pulses having a rise 
time characteristic of the mercury switch and an 
exponential decay of 2.5 m/xsec were used. It is in- 
teresting to note that when the two input pulses are 
of different amplitudes the response of the circuit is 
nearly independent of the amplitude of the larger 
pulse. From the curve of figure 2 corresponding to 
1.3-v output, it can be seen that the ratio of the 
amplitudes of a pair of equal coincident pulses to the 
amplitude of a single pulse producing the same re- 
sponse is about 1:25. The width of the curves of 
voltage versus signal delay is determined by the 
duration of the input. Higher pulses appear wider 
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Figuee 3. Response of the coincidence circuit of figure 1 and 
associated amplifier as a function of the time separation of a 
pair of equal amplitude pulses from the mercury-relay pulse 
generator. 

The dashed lines refer to operation of the circuit as part of a differential coinci- 
dence system (see text) . 

because of their finite rate of rise and exponential 
decay. 

In the experiments to be described, the photo- 
multiplier voltages were adjusted so that pulse 
amplitudes from 1 to 5 v were produced at the co- 
incidence circuit as a result of absorption in the 
scintillator of Compton electrons having energies 
between 200 and 1,000 kev. 

4. Coincidence System 

In the detection of gamma rays with organic 
scintillators, the energetic electrons produced in the 
Compton scattering of the gamma ray excite the 
phosphor, and the process of deexcitation results in 
the emission of light, which is detected by the photo- 
multiplier. The amount of light yield due to the 
interaction of a particular gamma ray depends on the 
energy of the scattered electron. Thus, by the nature 
of the Compton process, a monoenergetic gamma ray 
produces a continuum of pulse amplitudes at the out- 
put of the photomultiplier, with maximum pulse 
height depending on the gamma energy. 

As the response of the coincidence circuit depends 
on the amplitude of the input pulse, it is desirable to 
work with pulses above some minimum, amplitude. 
This can be accomplished without affecting the speed 
of response of the coincidence circuit by gating its 
output with the pulses from an amplitude discrim- 
inator into which the amplified pulses from, the last 
dynode of the photomultiplier are fed. Because the 
pulse distribution from the Compton electrons will 
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have a well-defined maximum, determined by the 
energy of the original gamma ray, one can study the 
response of the system to as narrow a band of pulse 
amplitudes as desired in the neighborhood of the 
maximum with the use of an integral pulse-height, 
selector. If one desires to work with pulse ampli- 
tudes below the maximum of the Compton distribu- 
tion, a differential pulse-height analyzer is required. 
For this study, integral pulse-height selectors were 
used with two sources, Na 22 with annihilation radia- 
tion coincidences of 0.510-Mev gamma rays and Co 60 
with coincidences between a 1.17- and a 1.33-Mev 
gamma ray. 

The resolution of the coincidence system as a whole 
is dependent on two factors. The first is the time 
distribution of pulses produced by the two photo- 
multiplier tubes when the original gamma rays are 
truly coincident. The second, factor is the duration 
of pulses from the photomultiplier tube because the 
coincidence circuit will respond to pulses that are 
delayed, relative to one another if the delay is less 
than the duration of the pulses. 

The usual procedure for studying the response of a 
coincidence system is to set up two scintillation 
counters, to excite the phosphors with coincident 
gamma rays such as those resulting from annihila- 
tion of positrons, and. to study the number of co- 
incidences recorded as a function of an artificial 
delay inserted in the system by varying the length of 
cable from one of the detectors to the coincidence 
circuit. If the pulses were square, of fixed duration, 
and. always coincident in time, the coincidence rate 
would, be constant as the cable length was varied 
over a range 4 corresponding to twice the duration of 
the pulse, and would be zero outside this range. 
On the other hand, if the pulse width were negligible, 
the curve of coincidence rate versus inserted delay 
would correspond, exactly to the time distribution 
of the pulses from, the detectors. 

The first step toward obtaining good resolution 
was the choice of the organic scintillating materials. 
One of the fastest known crystalline phosphors is 
trans-stilbene, which responds to a Compton electron 
by emission of an exponential pulse of light with a 
decay constant of about 6X10~ 9 sec. This was the 
phosphor that was chosen for the initial tests. 
Next, a series of tests was run on a number of RCA 
5819 photomultiplier tubes in order to select a few 
that were capable of producing pulses large enough 
to actuate the coincidence circuit without amplifi- 
cation. Some of the tubes tested, were capable of 
operation in the neighborhood of 2,000 v, and ex- 
hibited a maximum gain about 10 times that of the 
average tube. Two of those better tubes were 
chosen for use in the coincidence system. In this 
way, rise-time limitations introduced by amplifiers 
were eliminated. The duration of the pulse at the 
output of the photomultiplier tube can be reduced 
by using a shorted stub of coaxial cable to produce a 
delayed canceling pulse of opposite sign. The mini- 
mum, pulse duration obtainable by this procedure is 
approximately equal to the rise time of the original 
pulse. This rise time depends primarily on tin 4 
photomultiplier tbue, and appears to be 2 to 3X10 -9 



sec if the time constant of the output circuit is less 
than this value. 

As the best resolutions commonly attained with 
scintillation counters are of this order of magnit ude, it 
appears that the duration of the input pulses and 
fluctuations therein may often be limiting factors 
that mask the effects of the 4 time fluctuations in the 
pulses produced by coincident radiation events. 
A glance at figure 3 suggests a method to further 
reduce the effect of pulse duration upon the resolution 
of the s}^stem. Because of the width of the pulses, 
the rate of change of output voltage with delay is 
small in the region around zero delay. For this 
reason it is desirable to operate the circuit with such 
cable lengths that simultaneous pulses are delayed 
relative to one another by about 1 m./isec so that 
small fluctuations in the time of occurrence of pulses 
will produce a measurable change 4 in the output 
voltage. In order to have a symmetric system, it is 
necessary to use two fast coincidence circuits operat- 
ing on opposite sides of the maximum of the coin- 
cidence resolution curve, as indicated by the vertical 
dashed lines in figure 3. The outputs of these two 
circuits are fed to a slow coincidence 4 circuit. This 
complete 4 unit is then referred to as a differential 
coincidence 4 e'ire'uit, in analogy with a system pro- 
posed by Bay [2]. The 4 two individual fast coinci- 
dence circuits will be referred te> as channels of the 4 
differentia] e'oine*iele 4 ne'( 4 circuit. 

A coincidence is ree 4 ord.e 4 el only when there 4 are 4 
simultaneous output pulses from the 4 (we) fast coin- 
cidence circuits. The 4 cabling from the 4 dete 4 cte)rs is 
se) e'he)se 4 n that in the 4 event of a true 4 e'e)ine 4 ieh 4 neT the 4 
pulse from one 4 detector arrives before the 4 pulse 
from the other detector at channel 1 and after the 
pulse from the othe 4 r e!ete 4 e'te>r at e'hanne 4 l 2. 

As an illustration, it is useful to ce>nsiele 4 r the re- 
sponse of the circuit to square pulse 4 s of fixed ampli- 
tude. By making the differences in transmission 
time just less than the duration of the 4 pulse, all 
incidents of simultaneity of initial pulse 4 s result in a 
coincidence being recorded, whereas e 4 ve 4 n a delay 
that is less than the 4 duration of the 4 pulse 4 s results in 
a rejection by the differentia] circuit. The situation 
is shown schematically in figure 4. If the pulse from 
detector 2 is early, there will be no coincidence re- 
corded in channel 1, and if it is late 4 , there' will be no 
coincidence recorded in channel 2. Both channels 
will produce an output pulse only in case 4 of a true 
coincidence. 
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Figure 4. Principle of the differential coincidence technique as 
used in this investigation. 
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A block diagram of the differential coincidence 
circuit is shown in figure 5. The response of the 
system of coincident radiation is illustrated in figures 
6,7, and 8. The curves of figure 6 show the response 
of the individual channels to coincidences from co- 
balt-60. The displacement of the curves is to be 
attributed to the cabling arrangement between 
channels. The coincidence circuits are 100-percent 
efficient for pulses above the levels of the gating 
pulse-height selectors. The width of the curves is 
principally due to the width of the pulses from the 
photomuftipliers. 

Figure 7 shows the response when the two circuits 
are operated as a differential unit with the inserted 
delays such that there is no loss in efficiency. The 
fact that there is no loss in efficiency indicates that 
for zero delay the pulses overlap sufficiently at both 
channels to compensate for any time fluctuations 
that are introduced by the detectors. If the differ- 
ential cabling between channels is increased, there 
will be a decrease in the maximum coincidence rate 
when the time during which both pulses are present 
at the inputs of the coincidence circuit becomes less 
than the fluctuations in time of occurrence of the 
pulses. Figure 8 is a response curve for the system 
when this differential delay is sufficient to reduce 
the maximum coincidence rate to 10 percent of its 
previous value. The shape of this curve is deter- 
mined primarily by the time fluctuations of pulses 
from the detectors, and its width is a measure of the 
magnitude of these fluctuations. This conclusion is 
confirmed by the fact that a reduction in efficiency 
below about 20 percent results in no appreciable 
further narrowing of the response curve, indicating 
that the "channel width" in time is then much nar- 
rower than the distribution in time of the pulses 
arriving at the coincidence unit. A measure of the 
channel width is the quotient of the random coinci- 
dence rate by the product of the individual counting 
rates of the detectors, which, under the conditions 
of figure 8, is 2.8 X 10" 10 sec. 

Another factor that contributes to the width of the 
response curve is the variation in amplitude of the 
pulses at the input of the coincidence circuit. This 
effect has been minimized in this study by limiting 
the amplitude variation to about 25 percent with the 
gating pulse-height selectors. 
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Figure 5. Block diagram of the differential coincidence system 
with gating pulse height selectors. 
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Figure 6. Response of the individual channels of the coinci- 
dence system of figure 5 to pulses produced by coincident Co 60 
gamma rays on stilbene with 5819 photomultipliers. 

Different fixed lengths of cable interconnecting the inputs of the two circuits 
account for the displacement of the centroids of the two curves. The delay is 
introduced by varying the cable lengths from the photomultipliers. 



IUUU 
















































































































































































































































Id 




























































Z 










i 




















5 










/ 




























/ 




















-v 










/ 




















in 






























Id 
O 






























Z 
Id 
Q 






























O 

Z 

o 
° in 






































/ 






























1 






























I 




























5 




















































































i 








































I 












I 








I 










1 













-60 



-40 



-20 

DELAY X 10 ' 



20 40 

, SECONDS 



Figure 7. Response of the differential coincide nee system of 
figure 5, with cabling so arranged that the system responded to 
all coincident gamma rays producing pulses with amplitudes 
above the levels of the gating pulse-height selectors. 

5819 photomultipliers and stilbene crystals were used with Co 60 7-radiation. 
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Figure 8. Response of the differential coincidence system of 
figure 5, with cabling arranged so that the system responded to 
10 percent of the coincidences between gamma rays producing 
pulses above the level of the gating pulse-height selectors. 

The width of this curve is a measure of the time fluctuations in the scintillation 
counter. 
5819 photomultipliers and stilbeno crystals were used with Co 00 -y-radiation. 

5. Time Fluctuations in the Scintillation 
Counter 

5.1. Theoretical Analysis 

The differential coincidence circuit has been used 
to study the distribution in time of pulses from the 
photomultipliers under various operating conditions. 
Two sources of coincident radiation were utilized. 
Na 22 , a positron emitter, was used to obtain the co- 
incident 0.510-Mev gamma rays of annihilation. 
The number of unwanted coincidences between the 
1.3-Mev gamma rays also emitted by this source and 
the gamma rays of annihilation was made negligible 
by reducing the solid angles intercepted at the source 
by the detectors. Because the gamma rays of 2- 
quantum annihilation are emitted in exactly oppo- 
site directions, their coincidence rate depends only 
on a smaller of the two solid angles, whereas the rate 
of the coincidences involving the 1.3-Mev gamma 
rays diminishes as the product of the solid angles 
of the two detectors. A source of coincidences be- 
tween higher-energy gamma rays is Co 60 , which 
emits a beta particle followed by 1.17- and 1.33-Mev 
gamma rays in coincidence. 

Because the apparent time distribution of pulses 
from the detectors is determined by fluctuations in 
the rise times of the pulses, the factors that influence 
the rise time will be the factors that influence the 
observed distribution. These factors are fluctuations 



in photon emission and collection time and fluctua- 
tions in transit time of electrons in the photomulti- 
plier. Post and SchifT [3] have discussed the statis- 
tical limitations on the resolving time of a scint illation 
counter imposed by the fluctuations in photon arrival 
time, and Morton [4] has reported the results of 
a study in which account also has been taken of 
the transit-time fluctuations in the photomultiplier. 
Following the treatment of Post and Schiff, we con- 
sider a phosphor which, upon excitation at £=0, emits 
photons in an exponential decay process, so that 



><«=I' 



R\e~ xt dt, 



where f(t) is the average number of photoelectrons 
produced at the photocathode in the time interval 
0— t, R is the total number of photoelectrons pro- 
duced as a result of the excitation of the phosphor, 
and X is the decay constant of the phosphor. 



f(f)=R(X-e-*)siR\t f 



where £<C1/X. 



The number, /(/,), will be subject to statistical 
fluctuations described by the Poisson distribution, 
so that the probability of N electrons being produced 
in the interval 0— £ is 



Pn(0= 



IMfe- 



f(t) 



The probability that the Qth photoelectron is pro- 
duced between t and t+dt is 

W Q (t)=P Q -i(t)[df/dt]dt. 

The variance (v) in time of production of the Qth 
photoelectron is t 2 — (I) 2 . 

v= r t 2 W Q (t)dt~ ( °° tW Q (t)dtl 



«-«! 



-dt 



^Q/fix) 2 . 



Post and Schiff have considered the problem more 
generally, and obtain, for A£<1, 

<Q)= Wf [l+^^+liiglier-order terms in |]- 

To develop a model with which to compare the 
experimental results, Q was designated as the average 
number of photoelectrons that would be produced in 
a time equal to the average rise time of the output 
pulse. It is assumed that the variance in the time of 
response of the system to the pulse is given by the 
arithmetic mean of the variances associated indi- 
vidually with the first Q electrons. This assumption 
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Is reasonable, because the pulses are essentially 
differentiated by the use of the shorted stub, and the 
differential coincidence technique involves a coin- 
cidence between the leading edge of one derivative 
pulse and the trailing edge of the other. The average 
variance will be called Vi(Q). 

~^ Q) =uh[hM N+ i N+N2) }] 

1 V Q(Q+D | 2(Q+l)(Q+2) -\ 
~(i?X) 2 L 2Q ^ dR J 



Q 



2<2 2 



= 2(RX) 2 ' 3R(E\) 2 



where Q»l. 



The time of output of the signal will also depend 
on the fluctuations in electron transit time in the 
photomultiplier. Assuming, as did Morton, that in 
each stage the transit times can be described by a 
Gaussian distribution: 



,/ (t { -t ot y \ 
\ m f 



P(« )=,_£= exp 



where t oi is the average transit time in the itli stage, 
and Pi is a constant determined by the geometry and 
voltage in the iih stage. Let n t be the average 
number of electrons in the ith stage due to Q electrons 
at the photocathode. The variance in the average 
transit time of these electrons will be $/%. If <j is 
the average gain per stage, then n i = Qa (i ~ 1) . The 
average total transit time of the secondary electrons 
due to one photoelectron is 2t 0i , and the variance in 
this quantity, when Q photoelectrons are involved, is 



— ^^\ Pi ^r~\ Pi 



1^ $ 



For the average photomultiplier c~4, so that 
terms beyond the second can be reasonably neg- 
lected. Thus, the fluctuations in time of the output 
pulse of the photomultiplier are determined pri- 
marily by the transit-time dispersions in the regions 
between the cathode and the first dynode and be- 
tween the first and second dynodes. 



b(« 



Q 



+• 



In addition, there will be an apparent time fluctua- 
tion due to fluctuations in gain, primarily in the first 
stage. The apparent variance in time produced by 
the gain fluctuation will be denoted by v%=y 2 /Q, 
where 7 is a constant determined by the distribution 
of multiplication factors over a large number of 
processes. 
_To obtain the apparent variance in pulse time 
(F), we assume, as would be true for the convolution 
of Gaussian distribution, that the variance resulting 
from the several factors involved is the sum of the 



individual variances, so that 

^=^1 + ^2 + ^3 

Q , 2Q 2 



2(RX) 2 ' SR(R\) 2 



i(^7+" 2 ) 



a) 



where T is a constant characteristic of the photo- 
multiplier. 

Considering the case where Q/R is small, w^e can 
easily find the value of Q corresponding to minimum 
dispersion. 

* -LP o 



dQ=2(R\y Q* 
Q=^2R\T 

^2R\T . T 2 



V2T 



2(R\f [ J2R\T R* 



For a coincidence system using two photomulti- 
pliers, the variance of the response curve will be 
2v, and the width at half maximum of this curve 
of coincidence rate versus delay will be given by 



W=2A(2v) 



,-V/2 



(2) 



When QIR is small, we can utilize the expression 
for w mln to obtain 



W min ^4. 



\w 



(3) 



In practice we shall find that Q/R = 0.3 to 0.5, so 
that the value obtained from eq (3) for W m i n will 
represent a lower bound — perhaps 10 to 20 percent 
lower than the minimum obtainable width for the 
resolution curve. 

In order to compare the result displayed in eq (2) 
with experimental data obtained with a differential 
coincidence system, it is necessary to decide what 
factors determine Q for a system of this type. The 
shorted stub technique of pulse shaping will provide 
a pulse for the coincidence circuit where width is 
approximately equal to the rise time of the original 
current pulse at the anode of the photomultiplier. 
This rise time will be determined by the dispersion 
in the multiplier. We are interested in the number 
of photoelectrons (Q) that produce the secondary 
electron arriving at the anode during the rise of the 
pulse. As the rise time is determined by the dis- 
persion in the tube, the number (Q) will be propor- 
tional to the rate of production of photoelectrons, 
and we can write 

Q=KR\=K^2R\T. 
When Q/R is small, the dispersion is determined by 



v=K 



.,V2T 1 V2T = (K r2 +1) V27 7 
2R\ r K'2R\~~~ 2K' R\' 
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and the width at half maximum of the resolution 
curve is 



\ 2K' ) \R\) 



This expression exhibits a broad minimum in the 
neighborhood of K' = l, and if l/2<if / <2, the 
width of the resolution curve will be within 10 per- 
cent of the minimum obtainable with the combina- 
tion of tube and phosphor. 

It is of interest to consider in more detail the 
determination of the value of K f . The current 
pulse at the output of the photomultiplier can be 
described by 



■Hi'M^HM-*.}] 
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Where the zero of t and t m corresponds to the average 
time of arrival at the anode of the secondaries from 
a photoelectron emitted at the time of the inter- 
action in the phosphor, T measures the dispersion 
of the photomultiplier, and X is the decay constant 
of the phosphor. 

It can be shown in a straightforward manner that 
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We desire to obtain an expression for the rise 
time (/ r ) of this current pulse. At this time dg/dt=0. 
This condition yields the equation 
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(4) 



If X and T are known, this equation can easily be 
solved for t r by numerical procedures with the aid of 
Tables of Normal Probability Functions [5], which 
tabulates both functions. Best agreement with 
data presented in the next section is obtained by 
assuming a value of T=1.8X10 -9 sec in normal 
operation of the photomultiplier at 1,600 v. 

5.2. Comparison with Experiment 

Calculations have been made of the expected 
values of (t r ) for stilbene and terphenyl-toluene 
solution, assuming the above value (1.8X10 -9 sec) 
for Tand 1/X=6X10" 9 sec and 1/X=2.5X1CT 9 sec. 
respectively. For the two cases, 

t T =1.27T (for stilbene) 

t T = 0.867 > (for terphenyl-toluene). 



The maximum current occurs a time t r after the 
average arrival time of the secondaries produced by 
a photoelectron ejected at the time of interaction in 
the phosphor. The photoelectron charge, Q, con- 
tributing to the rise of the pulse is, on the average 



Q=R{ l-exp(-\/ r ) 



which yields 



Q=0.76\ r 2n\T (for stilbene) 
Q=0A5tJ2B\T (for terphenyl-toluene). 

In the case of stilbene, we neglect the term in Q 2 
in eq (1) to obtain from eq (3) a lower limit on the 
width at half maximum of the coincidence response 
curve. , We assume, in accordance with the work of 
Sangster [7], that stilbene gives 50 percent of the 
photoelectron yield of anthracene in a 5819 photo- 
multiplier. The yield of anthracene is half that of 
Nal measured by Hofstadter [8] to be about one 
photoelectron per kilo electron volt. For cobalt-60 
gamma rays, the average energy of the Compton 
electrons producing pulses accepted by the pulse 
height selectors is 900 kev. For stilbene we expect, 
therefore, 225 photoelectrons on the average. Using 
this value for R, we obtain a lower limit of 8.8 X 10 l0 
sec for the width at half maximum of the resolution 
curve. The expected value given bv formulas (1) 
and (2) is 9.7X10" 10 sec. The results of a measure- 
ment with the differential coincidence system are 
shown in figure 0, and confirm this value to within 
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Figure 9. Effect on the response of the differential coincidence 
system of figure 5 produced by a variation in the number of 
photons emitted by the phosphor. 

5819 photomultipliers and stilbene crystals were used. 
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10 percent, which is as good as can be expected in 
view of the uncertainties in the parameters involved. 
Figure 9 also shows the results of a measurement 
with annihilation radiation. In accordance with 
formula (2), we would expect a ratio of widths at 
half maximum of 



Wj 



Annihilation 



W c 



V^Annihilation/ 



-1.75 



because the average energy corresponding to the 
accepted pulses from the annihilation spectrum is 
300 kev. Using this ratio and the expected width 
of the resolution curve for cobalt-60 radiation, we 
find that the expected width of the low-energy curve 
is 17X10" 10 sec The observed width is 16X10" 10 
sec. 

In the case of the terphenyl-toluene scintillator, 
it is desirable to include all terms of eq (1) in esti- 
mating a lower limit for W, as Q/R is close to 0.5. 
The value of R is chosen in accordance with the work 
of Kallman and Furst [9] as 0.57 Rsmbene- For Co™ 
gammas, R=130 is used. The expression (1) for v 
can be minimized for our particular value of XT 7 
and, with (2), yields W 7 ml n=8.9X10- 10 sec. The ex- 
pected value of W is given by letting Q=0A5 ^2R\T 
in the expression for v. This gives W r C aic = 9.1 X 10~ 10 
sec. The experimental curve is shown as a dotted 
curve in figure 11, and yields W^bs^lOXlO -10 sec. 
For annihilation radiation, R= 130/3. The value of 
W calculated from eq (1) and (2) is 16X10" 10 sec. 
The experimental curve is shown in figure 10 along 
with that for stilbene, and gives W T O bs = 15X10~ 10 sec. 




-io o 

DELAY X 10" 

Figuee 10. Effect on the response of the differential coincidence 
system of figure 5 produced by a difference in the time distribu- 
tion of photons incident on the photomultiplier using annihila- 
tion radiation. 



In order to observe the effect of changing T, xhe 
voltage across the cathode to first dynode gap was 
increased by a factor of 4 and that from dynode 1 to 
dynode 2 by a factor of 2. From the geometries 
and voltages involved, it appears that £i~j3 2 . Ne- 
glecting 7, we have 



T\ 



*+? 



1 



because the quantity /3 is a measure of the mean- 
square spread in transit times and should be inversely 
proportional to the square root of the voltage. Thus 
the value of T 2 is (1.8X10- 9 )/V3J3 = 9.5X107 10 sec . 
Using this value for T 7 , we can calculate with the 
aid of eq (4) the rise time of the pulse on the anode 
of the tube. We find 

£ r =1.16r 2 , 
Q=0MR\T 2 . 

T 2 is sufficiently short to neglect the second-order 
term in calculating the lower limit on the width of 
the resolution curve. For cobalt-60 gamma rays on 
terphenyl-toluene we obtain H T min =5.4X10~ 10 sec. 
The calculated value of W, using eq (1) and (2), is 
6.1X10- 10 sec. 

The observed value of W obtained from the exper- 
imental curve of figure 11 is 6.2X10" 10 sec. 
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Figure 11. Effect on the response of the differential coincidence 
system of figure 5 produced by a change in the transit time in 
the photocathode-to-dynode region of the photo multiplier. 

5819 photomultipliers and terphenyl-toluene scintillators were used with Co 60 7- 
radiation. 
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Table 1. Width at half -maximum of curves of coincidence rate 
versus delay 



Scintillator 


Energy 


Voltage 


Tim in 


TVcalc 


U'obs 




kev 


V 


sec 


sec 


sec 


Stilbene 


900 


1,600 


8.8XIO-10 


9. 7X1O-J0 


9.0X10-1° 


Do 


290 


1, 600 


15.5 


17.0 


16.0 


Terphenyl-toluene. 


900 


1, 600 


8.9 


9.1 


10.0 


Do 


290 


1, 600 


15.5 


16.0 


15.0 


Do- 


900 


2,400 


5.4 


6.1 


6.2 



The results of these measurements are summarized 
in table 1. The agreement between calculated and 
observed widths of the resolution curves is excellent 
in view of the uncertainty in the values for E and X. 
These values were obtained from the literature. The 
values for R are considerably more in doubt , but the 
relative values should be within 10 to 15 percent. 
The value of the single parameter, T, was adjusted 
to give best agreement between the calculated and 
experimental values of W. An error in the absolute 
value of R would be reflected in this quantity. The 
value of 1.8X10 -9 sec obtained from this quantity 
does not appear to be unreasonable. Moreover, the 
consistency of the results suggests that the calcula- 
tion accurately represents the behavior of the coin- 
cidence system, and the conclusion may be drawn 
that in all cases observed the experimental resolu- 
tion was very close to the minimum obtainable with 
the combination of phosphor and phototube used. 



6. Multiplex Coincidence System 

The attainment of systems having better resolu- 
tion awaits the development of new phototubes 
capable of operating at very high voltages in the 
early stages and the discovery of faster phosphors, 
or requires the use of a different combination of 
crystal and multiplier. The latter alternative is 
illustrated by the following suggestion. If there are 
fluctuations in time of the pulses from the detectors 
for truly coincident events, the influence of these 
fluctuations on data indicating delayed radiation can 
be reduced by using two phot omul tip hers viewing 
each phosphor and requiring a delayed coincidence 
to be registered simultaneously in two independent 
coincidence circuits. Delayed radiation will produce 
a delayed coincidence in both channels, whereas a 
randomly occurring fluctuation is much less likely 
to have the same value in both sets of detectors at 
the same time. In fact, it has been experimentally 
verified that the curve of counts versus time for 
coincident radiation events falls off as the product of 
the curves for the individual coincidence circuits. 
The drawback on this procedure is that if two photo- 
multipliers are used, a reflector must be removed, and 
the amount of light collected by each photomultiplier 
is less than if only one tube is used. This decrease 
in light intensity corresponds to a decrease in R and 
has the same effect as a decrease in energy of tin 1 
original radiation. That is, the fluctuations in the 
response of the individual coincidence systems are 
increased. 
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Figure 12. Block diagram of multiplex coincidence system devised to minimize the effects of the time fluctuations in the scintillation 

counters. 
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A system based on this principle has been designed 
and constructed. It consists of two phosphors, each 
viewed by two photomultipliers. In this system, 
RCA 6199 photomultipliers were used because of 
their high gain and low noise. Four of the six 
tubes available were operated successfully at 2,000 v 
and produced larger pulses than the selected 5819 
photomultipliers used previously. The pulses from 
each phototube are fed to the inputs of two fast 
double-coincidence circuits, so that four of these 
circuits are used in the system. To minimize the 
effect of pulse width, the cabling is arranged so that, 
for zero delay and coincident radiation, each pulse 
arrives earlier than its counterpart at one coincidence 
circuit and later at the other. A time fluctuation in 
either direction will then eliminate the coincidence 
in 1 of the 2 circuits. The outputs of the fast 
coincidence circuits drive four inputs of a sixfold 
coincidence circuit. The other two inputs are 
driven by the outputs of two pulse-height selectors, 
each of which analyzes the summed pulses from the 
last dynodes of the photomultipliers viewing one of 
the phosphors. 

A block diagram of this system is shown in figure 
12. The system has been tested by means of the 
coincident gamma rays from Co 60 and from positron 
annihilation. The resolution curves are shown in 
figure 13. The resolution obtained with this system 
is about the same as that obtained with one difreren- 




SECONDS 



Figure 13. Response of the multiplex coincidence system of 
figure 12 to pulses produced by coincident gamma rays incident 
on terphenyl in toluene liquid scintillators. 

Co 60 gamma rays; , annihilation radiation. 6199 photomul- 



tipliers were used. 



tial coincidence system and 5819 photomultipliers 
operated at 2,400 v. From the standpoint of sta- 
bility, however, the new system has been much 
more satisfactory than the old one, primarily be- 
cause of a noticeable reduction in fatigue of the 
photomultipliers. 

7. Applications of Coincidence Techniques 

The mode of operation of high-speed coincidence 
circuitry is determined by the specific requirements 
of the problem at hand. If high efficiency is re- 
quired in the detection of a coincidence, a sacrifice in 
resolution must be made because of the inherent 
time fluctuations in the detectors. Such a situation 
is encountered, for example, in the absolute calibra- 
tion of radiation sources by coincidence methods. 
In this case, the system must operate with a con- 
stant known efficiency. An efficiency of 100 percent 
has been found to be the easiest to measure and hold 
constant, and a coincidence system has been de- 
veloped for the absolute calibration of cobalt-60 
sources by the method of gamma-gamma coinci- 
dences operating at this efficiency. For this applica- 
tion, the differential coincidence technique was not 
feasible because, with terphenyl-toluene scintillators, 
the maximum fluctuations in pulse time are of the 
order of magnitude of the pulse width. Because of 
the fluctuations in pulse time, it was necessary to 
accept a resolution r=1.5 + 10 -9 sec in order to 
achieve 100 percent efficiency. Under these condi- 
tions we are able to determine the disintegration 
rate of millicurie sources of cobalt-60 to an accuracy 
of 1 percent. 

If the problem is primarily one of measuring short 
time intervals between radiations, on the other 
hand, and less than 100 percent efficiency is tolerable, 
it is preferable to operate the coincidence system in 
such a way as to provide maximum resolution. In 
this application, it is possible to measure time inter- 
vals considerably shorter than the width of the reso- 
lution curve for the system. Bay [10] and Newton [11] 
have discussed some methods for analysis of the data 
obtained in such measurements. It can be shown 
that if the time delays are the result of the exponen- 
tial decay of an intermediate state excited by the 
first radiation, the curve of counts versus artificial 
delay is given by 



C(d)=AS^ (~) [exp { (j^j } ][exp { -a, 



-°cd}]X 






where A is a constant that depends on the geometry 
and resolution of the system, S is the source strength, 
ol is the decay constant of the intermediate state, and 
the prompt coincidence resolution curve of the 
system for radiations of the same energy is de- 
scribed by C(d) = C e K ' d \ The function contained 
in brackets has been tabulated in [5] so that it is 
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possible to calculate curves to be compared with the 
experimental curve and to determine the value of oc 
that gives a best fit. From a prompt coincidence 
resolution curve such as the solid curve of figure 13, 
it is possible to measure lifetimes of the order of 
10~ 10 sec to an accuracy of the order of 10~ 11 sec. 
These techniques have been applied to data obtained 
from the annihilation of positrons in metals and 
differences have been determined of 4X10 -11 sec 
between mean lives to a precision of 1X10 -11 sec. 
Details of this work will be published elsewhere. 

8. Conclusions 

The data and analysis presented here indicate that 
the circuitry and techniques described enable one 
to obtain coincidence resolutions very close to the 
minimum to be expected for the combinations of 
scintillator and photomultiplier presently available. 
The analysis of the detectors and circuitry takes 
account of the fluctuations in emission time of the 
photons from the phosphor and assumes Gaussian 
distributions for the interstage transit times in the 
photomultipliers. The time of occurrence of a pulse 
is considered to be the average of the arrival times 
at the anode of the electrons that contribute to the 
rise of the output pulse. This analysis results in 
the satisfactory prediction of the magnitude of the 
resolution and the effect on the resolution produced 
by a variation of tube voltage, phosphor decay time, 
or energy of exciting radiation. The data indicate 
that the root-mean-square deviation of the transit 
time in the cathode-to-dynode region of the 5819 
photomultiplier is 1.8X10 -9 sec when the voltage 
across the gap is 180 v. 



Fluctuations in photon-emission time and 
electron-transit time limit the time resolution obtain- 
able with a scintillation counting system, so that the 
use of circuitry having too short a resolution results 
in reduced efficiency for the detection of coincident 
radiation. For 1-Mcv gamma rays the limiting 
resolution for 100-percent efficiency is approximately 
1.5X10" 9 sec. Operation of the system at reduced 
efficiency enables one to measure time intervals of 
the order of 10~ 10 sec between radiations to an 
accuracy of 10 percent. 
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